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I. INTRODUCTION 

Searches for CP violation (CPV) in charm meson 
decays provide a probe of physics beyond the Stan- 
dard Model (SM). Singly Cabibbo-suppressed (SCS) de- 
cays can exhibit direct CP asymmetries due to inter- 
ference between tree-level transitions and |AC| = 1 
penguin-level transitions if there is both a strong and a 
weak phase difference between the two amplitudes. In 
the SM, the resulting asymmetries are suppressed by 
0(\V c bV u b/V cs V us \) ~ 10~ 3 , where Vj are elements of the 
Cabibbo-Kobayashi-Maskawa quark-mixing matrix 
A larger measured value of the CP asymmetry could be a 
consequence of the enhancement of penguin amplitudes 
in D meson decays due to final-state interactions 0, Q , 
or of new physics d, Q • 
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The LHCb and CDF Collaborations recently reported 
evidence for a non-zero CP asymmetry in the difference 
of the time-integrated D° — > tt + it~ and D° — > K + K~ 
decay rates @, 0] ■ Searches for CPV in other SCS de- 
cays with identical transitions c —¥ udd and c — > uss are 
relevant to an understanding of the origin of cpv @-Qil. 



We present here a study of the SCS decay D + — > 
K + K~it + [111 ], which is dominated by quasi-two-body 
decays with resonant intermediate states. This allows 
us to probe the Dalitz-plot substructure for asymme- 
tries in both the magnitudes and phases of the inter- 
mediate states. The results of this study include a mea- 
surement of the integrated CP asymmetry, the CP asym- 
metry in four regions of the Dalitz plot, a comparison of 
the binned D + and D~ Dalitz plots, a comparison of 
the Legendre polynomial moment distributions for the 
K + K~ and K~it + systems, and a comparison of param- 
eterized fits to the Dalitz plots. Previous measurements 
by the CLEO-c Collaboration found no evidence for CPV 
in specific two-body am plit udes or for the integrals over 
the entire phase-space [121 ]. The LHCb Collaboration 
also finds no evidence for CPV in a model-independent 
search [l3j . 
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II. THE BABAR DETECTOR AND DATA 
SAMPLE 

The analysis is based on a sample of electron-positron 
annihilation data collected at and just below the energy 
of the T(AS) resonance with the BABAR detector at the 
SLAC PEP-II collider, corresponding to an integrated 
luminosity of 476fb _1 . The BABAR detector is described 
in detail elsewhere 14]. The following is a brief sum- 
mary of the detector subsystems important to this anal- 
ysis. Charged-particle tracks are detected, and their mo- 
menta measured, by means of the combination of a 40- 
layer cylindrical drift chamber (DCH) and a five-layer 
silicon vertex tracker (SVT), both operating within a 
1.5-T solenoidal magnetic field. Information from a ring- 
imaging Cherenkov detector (DIRC) and specific energy- 
loss measurements (dE/dx) in the SVT and DCH are 
used to identify charged kaon and pion candidates. 

For various purposes described below we use samples of 
Monte Carlo (MC) simulated events generated using the 
JETSET [2l| program. These events are passed through 
a detector simulation based on the Geant4 toolkit |22j . 
Signal MC events refer to D + — > K + K~tt + decays gener- 
ated using JETSET as well as D+ A+A^tt+7 decays 
generated using JETSET in combination with the PHO- 
TOS [23] program. In all cases when we simulate particle 
decays we include EvtGen. [24[ 



III. EVENT SELECTION AND D+ K + K~n + 
RECONSTRUCTION 

The three-body D + — > K + K~tt + decay is recon- 
structed from events having at least three tracks with net 
charge +1. Two oppositely charged tracks must be con- 
sistent with the kaon hypothesis. Other charged tracks 
are assumed to be pions. To improve particle identifi- 
cation performance, there must be at least one photon 
in the DIRC associated with each track. Contamination 
from electrons is significantly reduced by means of dE/dx 
information from the DCH. Pion candidates must have 
transverse momentum px > 300 MeV/c. For lower px 
values, tracks are poorly reconstructed. Also, for lower 
Pt, differences in the nuclear cross sections for positively 
charged and negatively charged particles can lead to 
asymmetries. We form the invariant mass of K + K~tt + 
candidates and require it to lie within 1.82-1.92 GeV/c 2 . 
The three tracks must originate from a common vertex, 
and the vertex-constrained fit probability (P v tx) must be 
greater than 0.5%. The momentum in the CM frame 
(pcm) of the resulting D candidate must lie within the 
interval [2.4, 5.0] GeV/c. The lower limit on pqm re- 
duces background from B decays by preferentially se- 
lecting e + e~ — > cc events; this has traditionally been the 
way to reduce combinatoric background due to B decays. 
To remove background from misidentified D* + — > D°ir + 
decays, we require m(K + K~n + ) — m(K~iT + ) — to(7t + ) > 
15 MeV/c 2 , where the pion and kaon masses are set to the 



nominal values fl5j . Finally, for events with multiple 
candidates, the combination with the largest value of P v tx 
is selected. We perform a separate kinematic fit in which 
the mass is constrained to its nominal value . The 
result of the fit is used in the Dalitz plot and moments 
analyses described below. 

To aid in the discrimination between signal and back- 
ground events, we use the joint probability density func- 
tion (PDF) for L xy , the distance between the primary 
event vertex and the D meson decay vertex in the plane 
transverse to the beam direction, and pcm, to form a 
likelihood ratio: 

R = Ps{pcM)P s (Lx y ) 

C Ps{pCM)Ps{L xy ) + P h (p C M)Pb(L X y)' [> 

Since the two variables have little correlation we con- 
struct the two-dimensional PDF as simply the product 
of their one-dimensional PDFs; these one-dimensional 
PDFs for signal (P s ) and background (Pb) are esti- 
mated from data. The background PDFs are determined 
from events in the D + mass sidebands, while those for 
the signal are estimated from events in the D + signal 
region after background is subtracted using estimates 
from the sidebands. The signal region is defined by 
the m{K+K-TT + ) interval 1.86-1.88 GeV/c 2 , while the 
sideband regions are the 1.83-1.84 GeV/c 2 and 1.90-1.91 
GeV/c 2 intervals. The selection on R,£ is adjusted to 
maximize signal significance, and the resulting signal is 
fairly pure (see Fig. [3] in Sec. IVII) . 

The reconstruction efficiency for D + decays is deter- 
mined from a sample of MC events in which the decay is 
generated according to phase space (i.e., the Dalitz plot 
is uniformly populated). To parameterize the selection 
efficiency, we use the distribution of reconstructed events 
as a function of the cosine of the polar angle of the D 
meson in the CM frame [cos(#cm)] and the m 2 (A'~7r + ) 
versus m 2 (K + K~) Dalitz plot. The selection efficiency 
is determined as the ratio of N-R CCO /NQ cn in intervals of 
cos(#cm) and separately in intervals of the Dalitz plot, 
where -/Vr oco is the number of selected events in an inter- 
val and ./Veen is the number of events generated in the 
same interval. The binned Dalitz-plot efficiency is pa- 
rameterized with a feed-forward artificial neural network 
(ANN) [l6| consisting of two hidden layers with three and 
five nodes. Use of an ANN procedure allows us to ade- 
quately model the efficiency near the edges of the Dalitz 
plots. The ANN efficiency function is tested by creat- 
ing separate training and validation samples, which are 
satisfactorily fit by the ANN. 

IV. CORRECTIONS TO SIMULATED EVENTS 

In order to describe accurately the reconstruction ef- 
ficiency, we apply corrections to the reconstructed MC 
events to account for known differences between simu- 
lated events and data. The differences arise in the recon- 
struction asymmetry of charged-pion tracks, and in the 
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FIG. 1: (color online) Charged pion tracking efficiency asym- 
metry (defined in Eq. @) as a function of the pion momen- 
tum in the laboratory frame determined from the decays of r 
leptons. The horizontal error bars indicate the range of pion 
momentum [17l |. 

production model for charm mesons. Differences in kaon 
particle identification efficiency have a negligible asym- 
metry effect since the K + and K~ are common to D + 
and D~ decays. 

To correct the production model used in the simu- 
lation, we construct the ratio of the two-dimensional 
Pcm versus cos(#cm) PDFs between data and simula- 
tion and apply this ratio as a correction to the recon- 
structed MC events before calculating the efficiency. For 
this procedure the signal PDF for data is background- 
subtracted while the signal MC events are weighted by 
the Dalitz plot amplitude-squared determined from data 

(see Sec. Em}. 

To correct for differences in the reconstruction asym- 
metry of charged-pion tracks, we use a sample of e + e _ — ¥ 
t + t~ events in which one r decays leptonically via 
— > p^v^Vt while the other r decays hadronically 
via r T — > h T h T h ± v T . We tag events with a single iso- 
lated muon on one side of the event and reconstruct the 
hadronic r decay in the opposite hemisphere. We re- 
fer to this sample as the "Tau31" sample. We further 
require two of the three hadrons to have an invariant 
mass consistent with the rho mass to within 100 MeV/c 2 . 
Due to tracking inefficiencies, tau decays to three tracks 
are sometimes reconstructed with only two tracks. We 
use the two-dimensional distributions of cos n7r and pr w „ 
(with respect to the beam axis) of the rho-decay pions 
for 2-hadron and 3-hadron events to determine the pion 
inefficiency and asymmetry. We allow for a different effi- 
ciency for positive and negative tracks (e±) by introduc- 
ing the asymmetry a(pLab) as a function of pion labora- 
tory momentum (pLab): 



a(pLab) 



£+(PLab) - g-(PLab) 
£+(PLab) + £-G?Lab) 



(2) 



The results for a(j?Lab) are shown in Fig. [T] the av- 
erage value for < p Lah < 4 GeV/c is (0.10 ± 0.26)%, 



which is consistent with zero [17j . We use linear interpo- 
lation between data points, or extrapolation beyond the 
first and last data points, to obtain the ratio of track- 
efficiency asymmetries between data and MC as a func- 
tion of momentum. This ratio is then used to correct 
track efficiencies determined from signal MC. 



V. INTEGRATED CP ASYMMETRY AS A 
FUNCTION OF cos(0 C m) 

The production of D + (and D~) mesons from the 
e + e~ — > cc process is not symmetric in cos(#cm); this 
forward-backward (FB) asymmetry, coupled with the 
asymmetric acceptance of the detector, results in differ- 
ent yields for D + and D~ events. The FB asymmetry, 
to first order, arises from the interference of the separate 
annihilation processes involving a virtual photon and a 
Z° boson. We define the charge asymmetry A in a given 
interval of cos(#cm) by 



A , lQ \ \ N D +/e D + - N D -/e D - 

A(cos (0cm)) = ir T 1 — -r; 7 

N D + e D + + N D -/e D - 



(3) 



where N D ± and e£,± arc the yield and efficiency, respec- 
tively, in the given cos(6>cm) bin. We remove the FB 
asymmetry by averaging A over four intervals symmetric 
in cos(#cm), i.e., by evaluating 



A(cos(6» C m)) + A(- cos(6» CM )) 



(4) 



The interval boundaries in cos(#cm) are defined as [0, 
0.2, 0.4, 0.6, 1.0]. The D ± yields are determined from 
fits to the reconstructed K^K^ir mass distributions, 
as described in Sec. IVII This technique has been used in 
previous BABAR measurements in both three-body and 
two-body decays fl~8l - [20j | . The weighted average of val- 
ues obtained using Eq. @ is A CP = (0.37±0.30±0.15)%, 
where the uncertainties are statistical and systematic, re- 
spectively, with a probability of 21% that the asymme- 
tries are null in all four intervals (Fig. [2]) . 



VI. D+ MASS FIT 

The K + K~ir + mass distribution is fitted with a 
double- Gaussian function with a common mean and a 
linear background (Fig. [3]), plus a function describing ra- 
diative decays D + — > K + K~ir + 'y. The PDF for radiative 
decays is obtained from the reconstructed mass distri- 
bution of K + K ~7T + 7 events selected at the generator- 
level in our MC additionally convolved with a Gaus- 
sian of width 2.26 MeV/c 2 , and accounts for 1.5% of 
the signal. The fit to data gives a D + mass value of 
1869.70 ± 0.01 MeV/c 2 , where the uncertainty is sta- 
tistical only. The signal region is de fined to lie within 
±2(7d+ of the peak, where <jd+ = V h a \ + (1 /i)°2 
is 5.04 MeV/c 2 , and contains a total of 227874 events; 
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FIG. 2: (color online) CP asymmetry as a function of 
j cos(#cm)|- The solid line represents the central value of Acp 
and the dashed lines the ±1 standard deviation statistical 
uncertainty, determined from a \ 2 fit to a constant value. 

o\{a-i) is the standard deviation of the first (second) 
Gaussian component and /i = 0.63 is the fraction of 
the signal in the first Gaussian component. Separate 
fits to the K + K~tt + and K + K~tt~ distributions yield 
N D+ = 1 13037 ±469 and N D - = 110663±467 events, re- 
spectively. The ratio of efficiency-corrected yields (N/e) 



is R 



- N D+ /e 



p+ / e n+ — 



1.020±0.006. This ratio is used to ac- 



count for remaining asymmetries that arise from physics- 
or detector-related processes, such as an insufficiently ac- 
curate simulation of the FB asymmetry or a residual de- 
tector asymmetry. Also, it is a less accurate measure of 
the asymmetry when the efficiency varies significantly as 
a function of cos(#cm), as for our experiment. 



VII. MODEL-INDEPENDENT SEARCHES FOR 
CP VIOLATION IN THE DALITZ PLOTS 

Model-independent techniques to search for CP viola- 
tion in the Dalitz plots are presented in Ref. 18]. The 
techniques include a comparison of the moment distri- 
butions, and the asymmetry in the D + and D~ yields 
in various regions of the Dalitz plot. We scale the D~ 
yields by the factor R described in Sec. IVI1 By apply- 
ing this correction, we remove residual detector-induced 
asymmetries and decouple, as far as possible, the search 
for CPV in the Dalitz plot from the search for CPV in- 
tegrated over the phase space, which was described in 
Sec. El We measure the CP asymmetry in the four re- 
gions of the Dalitz plot labeled A, B, C, and D in Fig. [4] 
We report the fitted yields, average Dalitz plot efficien- 
cies, and CP asymmetries in Table Q] 

We pursue a second technique in search of CPV , 




m(K + K"7r±) GeV/c' 




1.92 



m(K + K"7t±) GeV/c 2 



FIG. 3: (color online) Combined reconstructed invariant mass 
distribution m(K + K~ 7r ) and projection of the fit result. 
The points show the data, the solid curve the fit model, and 
the dashed curve shows the background PDF. The signal re- 
gion is indicated by the dashed vertical lines, and the sideband 
regions by the solid vertical lines. The lower figure shows the 
fit on a logarithmic scale with the radiative component of the 
signal PDF shown separately as a smooth curve. 



by measuring normalized residuals A for the efficiency- 
corrected and background-subtracted D + and D~ Dalitz 
plots, where A is defined by 



n(D+) - Rn(D~) 
~ y/cr 2 (D+) + R 2 a 2 (D-)' 



(5) 



with n(D + ) and n{D~) the observed number of D + and 
D~ mesons in an interval of the Dalitz plot, where a(D + ) 
and a(D~) are the corresponding statistical uncertain- 
ties. The results for A are shown in Fig. [SJ Note that 
the intervals for Fig. [5] are adjusted so that each interval 
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TABLE I: Yields, efficiencies, and CP asymmetry in the regions of the Dalitz plot shown in Fig. [4] For the CP asymmetry, the 
first uncertainty is statistical and the second is systematic. 



Dalitz plot region 


N{D+) 


<D+)[%] 


N(D-) 


e(D-)[%] 


A CP [%] 


(A) Below AT* (892) u 


1882 ± 70 


7.00 


1859 ± 90 


6.97 


-0.7 ± 1.6 ± 1.7 


(B) K* (892)° 


36770 ± 251 


7.53 


36262 ± 257 


7.53 


-0.3 ±0.4 ±0.2 


(C) 0(1020) 


48856 ± 289 


8.57 


48009 ± 289 


8.54 


-0.3 ±0.3 ±0.5 


(D) Above K* (892)° 


and 0(1020) 25616 ± 244 


8.01 


24560 ± 242 


8.00 


1.1 ±0.5 ±0.3 
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FIG. 4: (color online) D — > K + K~ir Dalitz plot and fit projections assuming no CPV, with the regions used for model- 
independent comparisons indicated as boxes. The A/B boundary is at rriKn = 0.6 GeV 2 /c 4 , the B/C boundary at niKn = 
1.0 GeV 2 /c 4 , and the C/D boundary at rriKK = 1.3 GeV 2 /c 4 . In the fit projections, the data are represented by points with error 
bars and the fit results by the histograms. The normalized residuals below each projection, defined as (Abata — JVmc)/VjVmc, 
lie between ±5a. The horizontal lines correspond to ±3<r. 



contains approximately the same number of events. We calculate the quantity \ 2 /{ v ~ 1) = (J2i=i ^ 2 )li v ~ 1)> 
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Normalized Residuals 

FIG. 5: (color online) Normalized residuals of the D + and 
D~ Dalitz plots in equally populated intervals (top) and their 
distribution fitted with a Gaussian function (bottom). 

where v is the number of intervals in the Dalitz plot. We 
fit the distribution of normalized residuals to a Gaussian 
function, whose mean and root-mean-squared (RMS) de- 
viation values we find to be consistent with zero and one, 
respectively. We obtain \ 2 — 90.2 for 100 intervals with 
a Gaussian residual mean of 0.08 ± 0.15, RMS deviation 
of 1.11 ± 0.15, and a consistency at the 72% level that 
the Dalitz plots do not exhibit CP asymmetry. 

The Legendre polynomial moments of the cosine of 
the helicity angle of the decay products reflect the 
spin and mass of the intermediate resonant and non- 
resonant amplitudes, and the interference effects among 
them 25]. A comparison of these moments between 
the D + and D~ two-body mass distributions provides 
a model-independent method to search for CP viola- 
tion in the Dalitz plot, and to study its mass and spin 



structure. We define the helicity angle Oh for decays 
D + — > (r — > K + K~)n + via resonance r as the angle 
between the K + direction in the K + K~ rest frame and 
the prior direction of the K + K~ system in the D + rest 
frame. For decays D + — » (r — > K~n + )K + via resonance 
r, we define Oh as the angle between the K~ direction in 
the K~tt + system and the prior direction of the K~tt + 
system in the D + rest frame. 

The Legendre polynomial moment distribution for 
order / is defined as the efficiency-corrected and 
background-subtracted invariant two-body mass distri- 
bution m(K + K~) or m(K~ir + ), weighted by the spher- 
ical harmonic Y~, [cos(6»h)] = y/2l + 1/47tP;[cos(0h)], 
where Pi is the Legendre polynomial. We define 
the two-body invariant mass interval weig ht W 2 {1) = 

(Ej wf) s - J2 k w< lk B )/( e *)i where w ij( w ik) is the value 

of Yi for the j th (/e th ) event in the i th interval and (e^) 
is the average efficiency for the i th interval. The su- 
perscripts S and B refer to the signal and background 
components, respectively. The uncertainty on is 

<? {l) = v/E>F ) 2 +ZMk B ) 2 /^) 2 - To study dif- 
ferences between the D + and D~ amplitudes, we calcu- 
late the quantities X\ for I ranging from zero to seven in 
a two-body invariant mass interval, where 

, = (W l il) (D+)-RW^(D-)) 
^af{D+)+R^f{D-) 

We calculate the x 2 / n dof over 36 mass intervals in the 
K + K~ and K~Tt + moments using 

i h h 

where p\ l12 is the correlation coefficient between X lt and 
ph i 2 _ (X^X^)-(X^)(X^) (g) 

V{^ i)2 )-(^ (Zi) )V(^ 2)2 )-(^ 2) ) 2 ' 

and where the number of degrees of freedom (ndof) is 
given by the product of the number of mass intervals and 
the number of moments, minus one due to the constraint 
that the overall rates of D + and D~ mesons be equal. 
We find x 2 /ndof to be 1.10 and 1.09 for the K+R- and 
K~n + moments, respectively (for ndof = 287), which 
corresponds to a probability of 11% and 13%, again re- 
spectively, for the null hypothesis (no CPV). 

VIII. MODEL-DEPENDENT SEARCH FOR CP 
VIOLATION IN THE DALITZ PLOT 

The Dalitz plot amplitude A can be described by 
an isobar model, which is parameterized as a co- 
herent sum of amplitudes for a set of two-body in- 
termediate states r. Each amplitude has a com- 
plex coefficient, i.e., A r [m 2 (K + K~), m 2 (K~n + )] = 
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J2 r M r e l ^F r [m 2 (K+K-),m 2 (K-T: + )} [2jH2j|, where 
M r and r are real numbers, and the F r are dynamical 
functions describing the intermediate resonances. The 
complex coefficient may also be parameterized in Carte- 
sian form, x r — M. r cos(j) r and y r — M. r sm(f) r . We 
choose the K* (892)° as the reference amplitude in the 
CP-symmetric and CP-violating fits to the data, such 



that Mk* 



(892)0 



1 and 



'K* (892)° 



0. 



Using events from the sideband regions (defined in 
Fig. [3]) of the D + mass distribution, we model the 
CP conserving background, which is comprised of the 
A > *(892)° and 0(1020) resonance contributions and com- 
binatorial background. The combinatorial background 
outside the resonant regions has a smooth shape and 
is modeled with the non-parametric fc-nearest-neighbor 
density estimator [HI]. The AT* (892)° and 0(1020) re- 
gions are composed of the resonant structure and a lin- 
ear combinatorial background, which we parameterize as 
a function of the two-body mass and the cosine of the he- 
licity angle. The model consists of a Breit-Wigner (BW) 
PDF to describe the resonant line shape, and a first-order 
polynomial in mass to describe the combinatorial shape. 
These are further multiplied by a sum over low-order Leg- 
endre polynomials to model the angular dependence. 



TABLE II: Resonance mass and width values determined from 
the isobar model fit to the combined Dalitz-plot distribution. 

Resonance Mass (MeV/c 2 ) Width (MeV) 

A* (892)° 895.53 ± 0.17 44.90 ± 0.30 

0(1020) 1019.48 ± 0.01 4.37 ± 0.02 

a (1450) 1441.59 ± 3.77 268.58 ± 5.28 

A5(1430)° 1431.88 ± 5.89 293.62 ± 3.83 

K* (1680)° 1716.88 ± 21.03 319.28 ± 109.07 

/o(1370) 1221.59 ± 2.46 281.48 ± 6.6 

k(800) 798.35 ± 1.79 405.25 ± 5.05 



Assuming no CPV, we perform an unbinned 
maximum-likelihood fit to determine the relative frac- 
tions for the resonances contributing to the decay: 
AT*(892)°, K*(U3Q)°, 0(1020), a (1450), 0(1680), 
if 2 *(1430)°, AT*(1680)°, ^(1410)°, / 2 (1270), / (1370), 
/o(1500), ^(1525), k(800), / o (980), / o (1710), and anon- 
resonant (NR) constant amplitude over the entire Dalitz 
plot. We minimize the negative log likelihood (NLL) 
function 



-21n£ 

JV 

-2j> 

i=l 



JJ e M c(x 1 ,x 2 )S(xi,x 2 )dxidx 2 
B(xi,x 2 ) 



(1 -p(mi)) 



JJ B(xi,x 2 )dxidx 2 



(9) 



where N is the number of events. The reconstructed D + 
mass-dependent probability p(m) is defined as p(rrii) — 
S(m-)+B(m ) ' wnere S(m) and B(m) are the signal 



and background PDFs, whose parameters are deter- 
mined from the mass fit described in Sec. I VI I x\ — 
m 2 (K + K~) and x 2 = m 2 (K~Tr + ), S(xi,x 2 ) is the Dalitz 
plot amplitude-squared, £mc is the ANN efficiency, and 
B(x\,x 2 ) is the CP-symmetric background PDF. 

The mass and width values of several resonances, in- 
cluding the K*(892)° and 0(1020), are determined in the 
fit (Table HI)) . The /o(980) resonance is modeled with an 
effective BW parameterization: 



A 



1 



/o(980) 



m z a - m 2 - im T p KK 



(10) 



determined in the partial-wave analysis of Df — > 

K + K~tt + decays [30], where pkk = 2p/m with p the 
momentum of the K + in the K + K~ rest frame, mo = 
0.922 GeV/c 2 , and T = 0.24 GeV. The remaining reso- 
nances (defined as r — » AB) are modeled as relativistic 
BWs: 



RBW(M A b) 



F r Fo 



M 2 - M 2 AB 



iT AB M r 



(11) 



where Tab is a function of the mass Mab, the momen- 
tum pab of either daughter in the AB rest frame, the spin 
of the resonance, and the resonance width Tr. The form 
factors F r and Fd model the underlying quark structure 
of the parent particle of the intermediate resonances. Our 
model for the K~tt + S-wave term consists of the k(800), 
the ^q(1430)°, and a nonresonant amplitude. Different 
parameterizations for this term [3ll . [32j do not provide a 
better description of data. The resulting fit fractions are 
summarized in Table Mil We define a x 2 value as 



X 



E 



N MCl ) 2 



(12) 



where A^bins denotes 2209 intervals of variable size. The 
i th interval contains Ni events (around 100), and A*MCi 
denotes the integral of the Dalitz-plot model within the 
interval. We find x 2 / ndof = 1-21 for ndof = 2165. The 
distribution of the data in the Dalitz plot, the projections 
of the data and the model of the Dalitz plot variables, 
and the one-dimensional residuals of the data and the 
model, are shown in Fig. 2J 

To allow for the possibility of CPV in the decay, reso- 
nances with a fit fraction of at least 1% (see Table IIII[) are 
permitted to have different D + and D~ magnitudes and 
phase angles in the decay amplitudes (.4 or A). We per- 
form a simultaneous fit to the D + and D~ data, where 
we parameterize each resonance with four parameters: 
Mr,(j>r,rcp, and A0cp- The CP-violating parameters 



are tq p 



and A0cp = 



In the case 



\M r \ 2 + \M r \ 2 

of 5-wave resonances in the K + K~ system, which make 
only small contributions to the model, we use instead 
the Cartesian-form of the CP parameters, Ax and Ay, 
to parameterize the amplitudes and asymmetries. This 
choice of parameterization removes or eliminates techni- 
cal problems with the fit. For these resonances we there- 
fore introduce the parameters x r (D ± ) = x r ± Ax r /2 and 
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TABLE III: Fit fractions of the resonant and nonresonant am- 
plitudes in the isobar model fit to the data. The uncertainties 
are statistical. 



Resonance Fraction (%) 



K (892) 


oi i tr i n oa 

zl.lo ± U.zU 


<p{w20) 


oo /to _i_ n 1 o 


A o (1430) 


OCT OO _l_ O O A 

lb. 61 ± Z.z4 


INK 


6.38 ± 1.82 


«:(800) 


t no _l_ A *^o 

7.08 ± 0.6.5 


ao(145U) 


0.84 ± 0.69 


-f ( C\OC\\ 

/o(980) 


O A T _l_ A OA 

2.47 ± 0.30 


-f / 1 o ta \ 

fo (1370) 


1 1^7 _l_ A 1 

1.17 ± U.zl 


0(1680) 


0.82 ± 0.12 


i?r(1410) 


0.47 ± 0.37 


/o(1500) 


0.36 ± 0.08 


a 2 (1320) 


0.16 ± 0.03 


/a (1270) 


0.13 ± 0.03 


7f 2 *(1430) 


0.06 ± 0.02 


J?* (1680) 


0.05 ± 0.16 


/o(1710) 


0.04 ± 0.03 


/a (1525) 


0.02 ± 0.01 


Sum 


97.92 ± 3.09 



y r (D ± ) = y r ± Ay r /2. The masses and widths deter- 
mined in the initial fit (shown in Table [TTJ) are fixed, 
while the remaining parameters are determined in the 
fit. In Table llVl we report the CP asymmetries, i.e., ei- 
ther the polar- form pair {rep, A<pcp) or the Cartesian 
pair (Ax r , Ay r ). Figure [6] shows the difference between 
the Dalitz-plot projections of the D + and D~ decays, for 
both the data and the fit, where we weight the D~ events 
by the quantity R described in Sec. lVIl It is evident from 
the figure that both the charge asymmetry of the data 
and fit are consistent with zero and with each other. 



IX. SYSTEMATIC UNCERTAINTIES 

We consider the following sources of systematic un- 
certainty: the R,£ selection, corrections applied to the 
MC, binning of the data in cos(#cm), and the Dalitz plot 
model. 

To evaluate the uncertainty due to the R,£ selection, 
we vary the selection such that the yield varies by at 
least ±1 standard deviation, and assign a systematic un- 
certainty defined by the largest variation with respect to 
the nominal value of the CP asymmetry. 

The uncertainty due to corrections of the production 
model in the simulation (described in Sec. IIVI) is evalu- 
ated by randomly sampling the correction factors from 
a Gaussian distribution using their central values and 
uncertainties as the mean and sigma, respectively. The 
efficiency is then re-evaluated and the fit is re-performed 
floating the CP parameters while keeping other parame- 



ters fixed. This entire procedure is repeated 50 times. We 
take the RMS deviation of the 50 fit values of the CP pa- 
rameters to obtain the systematic uncertainty estimate. 
The uncertainty due to the tracking asymmetry correc- 
tion is evaluated by comparing the measurement with 
two different corrections, namely the "Tau31" correction 
and the correction used in our analysis of D + — > Kgir + 
decays [20| . The average tracking asymmetry in the lat- 
ter analysis is (0.23 ±0.05)%, which is consistent with the 
result presented in Sec. IIVI after accounting for the differ- 
ent momentum spectra. We take the difference between 
the CP asymmetry central values using the two different 
tracking asymmetry corrections as the systematic uncer- 
tainty. 

The integrated measurement results from binning the 
data in cos(#cm)- To evaluate the effect of the binning 
in cos(#cm) for the integrated CP measurement, we vary 
the number of intervals and the interval edges, and mea- 
sure the CP asymmetry as the average asymmetry from 
a single forward interval and a single backward interval. 
Systematic uncertainties are determined from the differ- 
ence between the nominal central value and the value de- 
termined from the alternative methods. We report these 
uncertainties for the integrated measurement in Table fVl 
These systematic uncertainties are combined in quadra- 
ture to obtain the final result. 

To determine the model-dependent uncertainty on the 
Dalitz-plot CPV parameters, we remove resonances with 
fit fractions less than 1%, one resonance at a time, and 
repeat the fit. We change the standard value of the ra- 
dius parameter in the Blatt-Weisskopf form factor [28| for 
the intermediate resonance decay vertex from 1.5 GeV -1 
to 1.0 GeV -1 . We take the maximum variation as the 
model-dependent systematic uncertainty. Systematic un- 
certainties for the Dalitz plot fit CPV parameters are 
listed in Table ED 

Finally, we study possible systematic effects on the 
binned Dalitz-plot results presented in Sec. I VIII The 
nominal probability for the null CPV hypothesis is 72% 
for 100 intervals, while it is 42%, 62%, and 73%, re- 
spectively, for 25, 49 and 144 intervals. In comparison, 
changing the R£ selection, as described above, changes 
the nominal probability to 81%. 



X. SUMMARY 

In summary, we do not find any evidence for CP viola- 
tion in the SCS decay D+ -> K+K~tt + . The integrated 
CP asymmetry obtained using Eq. (TJ| is (0.37 ± 0.30 ± 
0.15)%. We find that the asymmetries in four regions 
of the Dalitz plot are consistent with zero, as listed in 
Table HI In addition, the D + and D~ Dalitz plots are 
consistent with no CP asymmetry with a probability of 
72% according to the analysis of the normalized residuals 
for the D + and D~ Dalitz plot divided into 100 equally 
populated intervals. Finally, we find no evidence for CP 
asymmetry in decays through various intermediate states 
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TABLE IV: CP- violating parameters from the simultaneous Dalitz plot fit. The first uncertainties are statistical and the second 
are systematic. 



Resonance 


TCP (%) 




0(1020) 

^ *(1430)° 

NR 

k(800) 

a (1450)° 


u. ( r i j\ i.i ) i 

0.35t° .82 ±0.60 
-9.40t^ ± 4.42 

-14.301^57 ± 5 - 98 
2.00l4;gg ± 1-85 

5-07±6.54 ± 9 - 39 


7.43±1J| ± 2.35 
-6.llt3.24 ± 1.39 
— 2.56±^:?7 ± 8.91 

2.1012.45 ± 101 

4.00±|;96 ± 3 - 83 


/o(980) 
/o(1370) 


Ax 

-0.199lo:uo ±0.084 
0.019l£o4s ± 0.022 


Ay 

-0.231±Jio5 ± 0.079 
— 0.0045±g;g|| ± 0.016 
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FIG. 6: (color online) The difference beween the D + and D Dalitz plot projections of data (points) and of the fit (cyan band). 
The width of the band represents the ±1 standard deviation statistical uncertainty expected for the size of our data sample. 



TABLE V: Systematic uncertainties for the integrated CP 
asymmetry. 



Average cos(#) asymmetry 


AAcp[%\ 


Event selection 


0.07 


Single forward and backward bin 


0.01 


cos(#cm) binning 


0.04 


Track asymmetry correction 


0.12 



from a study of the two-body mass distributions, as seen 
in Fig. [51 and from a parameterization of the Dalitz plot 
for which the CP asymmetries in amplitudes are listed in 
Table El 
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